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Abstract 

A general formalism for the evaluation of time reversal odd parton distributions is applied here 
to calculate the Boer-Mulders function. The same formalism when applied to evaluate the Sivers 
function led to results which fulfill the Burkardt sum rule quite well. The calculation here has 
been performed for two different models of proton structure: a constituent quark model and the 
MIT bag model. In the latter case, important differences are found with respect to a previous 
evaluation in the same framework, a feature already encountered in the calculation of the Sivers 
function. The results obtained are consistent with the present wisdom, i.e., the contributions for 
the u and d flavors turn out to have the same sign, following the pattern suggested analyzing the 
model independent features of the impact parameter dependent generalized parton distributions. 
It is therefore confirmed that the present approach is suitable for the analysis of time reversal odd 
distribution functions. A critical comparison between the outcomes of the two models, as well as 
between the results of the calculations for the Sivers and Boer-Mulders functions, is also carried 
out. 
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I. INTRODUCTION 



The study of the transverse polarization of quarks in the nucleon, one of its less known 
features (for a review, see, e.g., Ref. [l| ) , is progressing fast, since the expected new data from 
ongoing experiments are motivating a strong theoretical activity (for recent developments, 
see Ref. |2j). The present work aims to contribute to the understanding of the transverse 
polarization of quarks by evaluating, using different models of the proton structure, the 



Boer-Mulders function 3j. We will use a theoretical scenario 



the calculation of the Sivers function 



sign, the Burkardt sum rule and the general trend of the parameterizations extracted 
from the available data 8j. 

The Boer-Mulders function h± describes the number density of transversely polarized 
quarks in an unpolarized target. As the Sivers function, fyj,, describing the number density 
of unpolarized quarks in a transversely polarized target, the Boer-Mulders function is a 
Transverse Momentum Dependent (TMD) PD and it is a time-reversal odd object (T-odd) 
For this reason, for several years, it was believed to vanish due to time reversal invariance. 



recently developed for 



6] , which has reproduced its main features, such as the 



However, this argument was inva 
calculation in a spectator model 



idated, initially in the case of the Sivers function, by a 
9], following the observation of the existence of leading- 

n 

twist Final State Interactions (FSI) [10] . The current wisdom is that non- vanishing T-odd 
unctions are generated by the gauge link in the definition of TMD parton distributions 
11 . [12I whose contribution does not vanish in the light-cone gauge, as it happens for 

the standard PD functions. Those T-odd functions can be intuitively related to impact 

parameter dependent (IPD) parton distributions 

link, it is formally questionable. For example, in the case of the Sivers function, it has been 



[14], Q. However, because of the gauge 



shown that a model independent relation between this quantity and the corresponding GPD 
is still to be found [if]]. 

The Boer-Mulders function is being measured through Semi-inclusive deep inelastic scat- 



tering (sidis) 0, Q, Q (see Ref. 

Drell-Yan (DY) process in hadronic colliders 



20] for an analysis of these data) and through the 



21 



22 



23] . In particular, in Ref. 24|, the 



Boer-Mulders function has been recently extracted from the data of the unpolarized p + D 
Drell-Yan processes measured by the E866/NuSea Collaboration at FNAL |22|. However, 
the extraction of h± is very difficult, as it always involves another chirally-odd distribution 
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function. As a matter of fact, in SIDIS, the Collins fragmentation function is required 25], 
and in DY the observable quantity is a convolution of two belonging to the two colliding 



nucleons 



26|. Presently, the lack of accurate data affects the quality of the extraction. 



The present experimental scenario motivates therefore the formulation of theoretical esti- 
mates. In principle one should perform a calculation in QCD; however, this is presently not 
possible. Lacking this possibility, it becomes relevant to perform model calculations of the 
Boer-Mulders function. A few estimates exist. In a quark-diquark model with axial vector 
diquarks, it was originally found that has a different sign for the u and d flavors 27] . 



2s 



who, following a procedure 



These findings have been corrected by the Authors of Ref. 
established in previous papers 29| , demonstrated that the sign of for the u and d flavors 
turns out to be the same, assuming the spectator diquark to be either scalar or axial vector. 
This feature is in qualitative agreement with the pattern predicted by quark helicity-flip 
IPD GPDs in models 



30 



3jj and in lattice simulations [32] . It has also been found in a few 



other model calculations: in the MIT bag model, in its simplest version 



331 ] : in a large N c 



analysis of TMDs 34|; in a phenomenological parameterization based on the quark-diquark 



picture 



35]. 



4( and 



used so far for the evaluation 



In here, the recently proposed formalism of Refs. 
of the Sivers function in a Constituent Quark Model (CQM) and in the MIT bag model, 
respectively, will be extended to calculate hj- for the valence quarks. In the case of the 
Sivers function, within both models, this approach has proven to be able to reproduce the 
main features of /fy. Similar expectations motivate the present analysis of h^. The MIT 
bag model calculation presented here has an additional purpose, namely completeness. As 
a matter of fact, in Ref. 33|, the general framework for the calculation of T-odd TMDs 



has been nicely set up, but an important contribution has been nevertheless disregarded. In 
Ref. , we have reincorporated this contribution into the calculation of the Sivers function. 



We found that the results of Ref. [33j were incomplete and, as a consequence, they did not 
fulfill the Burkardt Sum Rule. Once this missing contribution is properly included in the 
calculation, the Burkardt Sum Rule turns out to be fulfilled to a large extent [5]. The same 



problem affects the evaluation of hj; in Ref. 33j , which is therefore retaken here, along the 
lines of Ref. fl. 

The paper is structured as follows. In the second section, the main quantities of interest 
are introduced and the formalism for the calculation of in a CQM and in the MIT bag 
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model is developed (some technical details are given in the Appendix). The numerical results 
of the calculations are presented and discussed in the following section. In the last section, 
we draw conclusions from our study. All throughout the paper, the formalism and the results 
are critically compared with those previously obtained for the Sivers function f^. 



II. FORMALISM 



The Boer-Mulders function, h t (x, kx) 3 



the quantity of interest here, is formally 



defined, according to the Trento convention 36j, for the quark of flavor Q, through the 
following expression 1 : 



2k x J (2vr) : 



x \ E (^I^Q(r,fr)4 T (oo,r)7V75A(oo,0)^Q(0,0)|P^) + h.c. , 

S z =-l,l 

(1) 



where S is the spin of the target hadron, the normalization of the covariant spin vector is 
S 2 = — 1, M is the target mass, "0s(£) is the quark field and the gauge link is 

%(oo,r) = Vexp(-ig J A+(r,-,i T )dr,-j , (2) 



where g is the strong coupling constant. One should notice that this definition for the gauge 
link holds in covariant (non singular) gauges, and in SIDISprocesses, since the definition of 
T-odd TMDs is process dependent. As observed in Ref. j^l for the first time, and later in 



37 



38] using factorization theorems, the gauge link contains a scaling contribution which 
makes the T-odd TMDs non vanishing in the Bjorken limit. 

For completeness, we recall now the definition of the Sivers function, f^r(x, fcy). Taking 



1 An equivalent expression is obtained by performing the following changes in Eq. (TT]): k x — > —k y and 



7 2 — > 7 1 . Here and in the following, = (ao ± 03) /v2 and fcr = \kr\- 
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the proton polarized along the y axis, one has 



-i(xp+£ -k T -£ T ) 



2k x J (2vr) 



l - S y (PS y \4j Q (r,&)Ci T (oo,C)l + C (oo,0)M0,0)\PS y ) +h.c. . 

Sy = -l,l 

(3) 



The difference between the two quantities is clearly seen and physically transparent in Eqs. 
(CQ) and ([3]). The function (x, kx) counts the transversely polarized quarks (as given by 
the Dirac Operator 7 + 7 2 75 in Eq. (CQ)) in an unpolarized proton (as given by the average 
over the proton helicity S z in Eq. (TjQ)), while /f^ S (x, kj) counts the unpolarized quarks (as 
given by the Dirac Operator 7 + in Eq. ([3])) in a transversly polarized proton (as given by 
the explicit transverse component S y in Eq. ([3])). If there were no gauge links, the two 
T-odd functions f^(x, kx) and h^ Q (x, kj>) would be identically zero. Expanding the gauge 
link to the first non-trivial order, i.e. the next to leading one, the Boer-Mulders function 
h± (x, kr), Eq. (CQ), becomes 

hl {X ' kT) ~ 2k x J ( 27 r)3 6 



X 2 

s*=-i,i 



/CO 
A+(r / -,f T )^-^7 + 7 2 75^(0,0)|P5 2 )+h.c. . (4) 

The formalism will be presented now for the evaluation of h^(x, kx) in the MIT bag model 
and in the CQM. Let us start from the former. At the beginning, our procedure will follow 



step by step the one nicely set up in Ref. 33| for the evaluation of T-odd TMDs in the MIT 



39| |) are given 



bag model. The properly normalized fields for the quark in the bag (see also 
in terms of the quark wave function in momentum space, which read 

<p n {%) = i V^NRl[ tom * m | , (5) 

\a- fcti(|fc|) Xn 



with the normalization factor N 

N = 



3 \ 1/2 



2RI(uj- l)sin 2 o; 



5 



where uo = 2.04 for the lowest mode and Ro = 4u/M is the bag radius. The functions U{k) 
are defined as 



U(k) 



u 2 duji(ukRo)ji(uuj) . 



(6) 



Following Ref. 33j, we fix the other ingredients of Eq. (j4]), having in mind Fig. 1 (and its 
h.c). In particular, the gluon propagator is treated in a perturbative way and it is assumed 
that it is not modified in the bag medium. This leads, in the Landau gauge, to 

g 2 M E P f d 3 k 3 d 3 ki d 4 q 



hf Q (x,k T ) 



k x 

x 27r<5(g ) 



(2vr) 3 (2tt) 3 (2tt) 4 
1 



5 (kf - xP + + q + ) 5 {2) (k 1T -k T + q T \ 



ie q z + it 



^EE E WM( PS *\4rnArn/L bl ernJPS Z ) 
Sz=— 1,1 P mi,m2,m3,i7i4 

x fLfii) 7°7 + 7 2 7s fefiife) 7°7 + ^m 4 (h - 0) + h.c. . 



(7) 



The last expression corresponds to the definition of the Boer-Mulders function in the MIT 
bag model, given in Ref. 



Substituting in Eq. (J7|) the identity 

1 



1 



g + — ie q + + ie 
and performing the trivial integrations, one gets 



i 2tt 8(q- 



d 3 k? 



Sz=— 1,1 /3 mi,m2,m3,m4 

x VmiC*- *r)7 7 + 7 2 75^m 2 (A;) 



7°7 + ^3 - <?t) • (8) 



(27r) 3 ^ m3 

The last line here depends on the four spin indices, whose combinations are dictated by 
the spin-flavor-color matrix elements given in the second line. Let us see now where the 
differences between the present calculation and that of Ref. l33| arise. Following the notation 

n 

of Ref. [5|] , we write 



flm $ - <?t)7°7 + 7 2 75 fm 2 (k) = Inn $, qr) S mi rn. 2 + J mi (k, q T ) 6 mi , 



-m,2 ) 



(9) 
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/d?k -> -> 

J^T^ < fL 3 (h)l°l +< fm4(h - Qt) = F m3 (q T )S m3m4 + H m3 (q T )5 m3 _ m4 . (10) 

In these expressions, all the possible helicity combinations in the initial and final states are 
emphasized. To our understanding, there are no physical reasons to rule out none of the 
above possible contributions. It is easily realized that both the terms with no-helicity flip 
of the quark 3, i.e. F m3 , as well as the one allowing the helicity-flip of this quark, i.e. H m3 , 
are non-vanishing under the integration over d 3 k^. This result is in contrast with the result 



presented in Ref. 33j, where only a term similar to F m3 is considered. As we explained in 
Ref . j^] , this result only applies if the integral is performed taking cfr along the z direction. 
However, this is incorrect because, in any DIS process, the operator structure is determined 
by the direction of the virtual photon. Here, the operator 7+ results from having taken 
the photon along the z axis. Therefore we no longer have the freedom to choose z as the 
direction of the exchanged gluon, which must then lie in the (x, y) plane. Besides, one 
can check that the integral Eq. ( flOl) does depend on the direction of qr- Moreover, if the 



findings of Ref. 33| were correct, it would mean that no helicity flip for the interacting quark 
could occur, a restriction which does not have any physical motivation. Thus the present 
calculation differs from the previous one in that we take into consideration both terms of 
Eq. (TTO]) . By the same argument, the expression tp^^k — qT)l°l + l 2 l^rn 2 {k) m Eq. (jSJ) also 
contains both helicity-flip and non-flip terms. 

Using Eqs. (J9J) and (fTUl) in Eq. (jBJ), the Boer- Mulders function reads 

2 ME P f d 2 q T 1 



hi Q (x,k T ) = -2W 2 -^ J 



(2tt)V 

S z =— 1,1 /3 mi,m2,m3,m4 

{k,q T )5 

1711,-1112 }{^ 3 7773,-7714 f • 

(11) 

After the evaluation of the spin-flavor-color matrix elements, 

performed assuming SU(6) symmetry, for Q = u,d and S z = 1(— 1) (these coefficients are 
listed in the Appendix), and after a straightforward calculation, one is left with the final 
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expression 



,-Lu(d), , x 2 ME P f d 2 q T C 

K '{x,k T ) = -g 



k x J (2tt) 2 q 2 

x {6(3) Ji F ± + 1(2) J 2 F 2 + 1(2) J\ H ± + 1(2) J 2 if 2 } , (13) 

with the eight functions Ii )2 , -Pi,2, #i,2> ^1,2 and the c-number C given in the Appendix. The 
two first terms on the r.h.s. of ffTB"]) are the contributions involving quarks which do not flip 
helicity. On the other hand, the last two terms on the r.h.s. of the same expression represent 
the contribution due to the double helicity flip of the quarks. This contribution will turn 
out to be non negligible. We reiterate that the result Eq. (Tl3l) is different from that of the 



previous calculation in the bag model [33j, which includes only a F m3 -like term. 

Through the coefficients Eq. ( 1121) it is possible to reconstruct what happens at the level 
of the quark helicity in a perfectly transparent way. The role of each contribution can 
be followed and evaluated. The dominant contributions are the non-flipping I\Fi and the 
double-flipping J±Hi ones. Due to the spin-flavor-color coefficients, i.e., due to the SU(6) 
symmetry assumption, the non-flipping term is bigger than the double-flipping contribution. 
The improvement of the Boer-Mulders in the MIT bag model with respect to the previous 



calculation [33|] quantitatively consists in the addition of the double-flipping contribution. 
In effect, the other two terms, I2F2 an d J2H2, are governed by the product of the two lower 
components of the bag wave function (cf. Eqs. (A8) and (A14) in the Appendix), which 
encodes the most relativistic contribution arising in the MIT bag model. They turn out to 
be a few orders of magnitude smaller than the dominant ones, arising from the interference 
between the upper and lower parts of the bag wave function. We will see now that this 
happens also if a proper non relativistic reduction (NR) of the gauge link, suitable for CQM 
calculations, is performed. 

The calculation scheme for T-odd functions in a CQM has been completely set up in 
Ref. [4], where it has been applied to the Sivers function. Since the formalism required here 
for the Boer-Mulders function is the same used in our previous paper to evaluate /f L2 (x, kx), 
we refer to that paper to obtain a workable formula for (x, kx). As it has been discussed 
after Eqs. (1) and (3), with respect to the fi(x,kx) calculation, one has to change the 
polarization of the proton and the Dirac structure of the quark operator. Starting from the 
definition, Eq. ([1]), it is found therefore that, in a CQM, the equivalent of the expression 
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Eq. (jSj), previously obtained in the MIT bag model, is: 

r2 r Ai, 



M C da 

hi Q (x, k T ) = -2ig 2 — / dhdh-^Sikt - xP + )5(k 1T + q T - k T )M Q , (14) 

k x J {^) 



where 



— ^ ] ^(5 z ,r,mi,m2,i„,m3,"i4,ij,C3,c 4 ) 2 ] 

J r 3,mi,C3,Tri2,C4,m3,i,m4,j S z =— 1,1 

x \&J gjs Ui{mi,i, Q}; fc 3 {^3, c 3 , .F 3 }; ~h ~ hjn^ 
x W^V h (k x ,k 3 ,q) 

x Ui + ?,{m 2 ,i, Q}; 4-9, {m 4 ,c 4 , ^3}; -£ 3 -£i»k) > ( 15 ) 

where, Zj = {mj,Q, Q} represents the set of helicity, color and flavor quantum numbers, 
respectively, and the vector is an intrinsic proton state. Given the diagram shown in 
Fig. 1, the interaction V h (ki,ks,q) reads 

V h (kx,k 3 ,q) = -jU mi (X) 7 + 7 2 75^m 2 + <f) ^m 3 (4) 7 + ^m 4 (X-<f) • (16) 

f r s z can be factorized into a completely antisymmetric color wave function, T, and a sym- 
metric spin-flavor-momentum state, $ s /, as follows: 

*rs 2 = ^{7713,^3}; -fei-^K,^}) r^c^Cn) . (17) 

The matrix element of the color operator in Eq. (1151) can be therefore immediately evaluated: 

C 3 , C n )T^ iC4 T(j, C4, Cn) = -\, (18) 

C3,C4,iJ 

which is the well-known result for the exchange of one gluon between quarks in a color 
singlet 3-quark state 40j. Besides, as a consequence of the symmetry of the state $ s /, one 
can perform the calculation assuming that the interacting quark is the one labeled "1" , so 
that, after the evaluation of the summation on the flavors Ai a can be written, for the u 
and d flavors, as follows: 

\ S Q — _1 1 m, rr>„ »no ™ , 



X 



Sz=— 1,1 mi,m2,ni3,m4 

1± 73(1) 



2 



x $ s/j5z (fa + <f, m 2 ; fc 3 - q, -fa - k 3 , m^j , (19) 



where the helicity m n of the spectator quark is determined by mi, m% and S z due to angular 
momentum conservation. Some remarks concerning the interaction, Eq. ffl6l) . are in order. 
In general, the interaction can be separated according to the possible helicity combinations 
of the interacting quarks in the initial and final state. Besides, since the wave functions to 
be used are NR, a NR reduction of the interaction has to be obtained. This NR reduction 
is given by the V^ R (ki,ks,q) expression in Eq. ([191) . In order to obtain it, we follow the 
procedure developed for the Sivers function in Ref . [4 1 , using therefore the definitions of free 
four-spinors in Eq. (fl6j) and performing a NR expansion as it is commonly done in nuclear 
physics. At O (k 2 /M 2 ) one obtains 

C (k u k 3 , g) = VZ~ fHp (4 h, q) + V N d ° uble - fHp (h, h, q) , (20) 

i.e., the relevant processes for the evaluation of h± are the ones where either the interacting 
and active quarks do not flip their helicities, or they flip them both. The expressions for 
the p and V^™ le ~ potentials are given in the Appendix. As explained in Ref. 4] 

in the case of the Sivers function, we reiterate that, in this approach, it is the interference 
of the upper and lower components of the four-spinors of the free quark states which leads 
to a no n- vanishing h^. Effectively, these interference terms in the interaction are the ones 
that, in the MIT bag model previously described, arise due to the wave function. Terms of 
higher order would correspond to those arising, in the MIT bag model calculation, from the 
product of the lower components of the quark wave function, a contribution which has been 
found to be negligible. The results of the relativistic calculation justifies therefore the NR 
reduction of the interaction in the CQM calculation. 

Eq. ( Tl4l) . with Ai u ^ given by Eq. (fl9l) . provides a suitable formula to evaluate 
using the proton state in momentum space, described in a CQM. We will here restrict 
our calculation to a Harmonic Oscillator potential model with pure SU(6) symmetry for 
the proton. The choice of a SU(6) wave function is motivated to render the comparison 
with the previous calculation, performed in the MIT bag model with SU(6) symmetry, more 
plausible. It is worth recalling that, as it has been stressed in Ref. 4J, the breaking of SU(6) 



in the CQM of Isgur and Karl 



4l| does not change essentially the results of the calculation. 



This has been obtained also in the evaluation of standard PDs or GPDs in quark models 



see, i.e. 42|): the essential features of a CQM calculation are obtained in a pure SU(6) 



framework. The formal expression of the proton state can be given in terms of the sets of 
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conjugated intrinsic coordinates given in the Appendix. It reads 

e -(kl + ki)/o* 



\ 2 S\/2)S 



■\xh 



(21) 



7r 3 / 2 a 3 

where the spectroscopic notation \ 2S+1 Xj) t with t = A,M,S being the symmetry type, has 
been used, and \x)s is the standard SU(6) vector describing the spin-flavor structure of the 
proton. The parameter a 2 = moj of the H.O. potential is fixed to the value 1.35 fm~ 2 , in 
order to reproduce the slope of the proton charge form factor at zero momentum transfer 



40] 



We have now all the necessary ingredients to write the final expression for h x 
reads 



u(d) 



which 



hf Hd \x,k T ) 



g 



,M 2 



3/2 



2tt 3 / 2 



1 / d 2 q T fl V2k° x 1 

! / 2 a 3 J (2vr) 2 V 2 \k° x - k z x \ q 2 6 " 



q.k; 



K 1 3 V 3 J I 48m 2 



(22) 



with k° x 



^Jm 2 + /c 2 , and 



k\ = \l ^(q~ k) 



\m 2 + k\ 



XT 



3x 2 P+ 2 



2 V / 3P H 



x 



kf + k 



XT 



(23) 



The explicit expressions for the functions A, B are given in the Appendix. 



III. RESULTS AND DISCUSSION 



To evaluate numerically Eqs. (1131) and (I2"2"j) . the strong coupling constant g, and therefore 
a s (Q 2 ), has to be fixed. Here, the model-independent prescription introduced for calculations 
of PDs in quark models (see, i.e., Ref. 42]) will be used. It consists in fixing the momentum 
scale of the model, the so-called hadronic scale /Xq, according to the amount of momentum 
carried by the valence quarks. In the present approach, for both the MIT bag model and 
the CQM, only valence quarks contribute. The assumption that all the gluons and sea pairs 
in the proton are produced perturbatively, according to NLO evolution equations, yields 
$ ~ 0.1 GeV 2 , if Ag£g ~ 0.24 GeV. Therefore, one has a s (/x§)/(47r) ~ 0.13 |42[ . This is 
obtained imposing that ~ 55% of the momentum is carried by the valence quarks at a scale 
of 0.34 GeV 2 , as in typical low-energy parameterizations (see, i.e., Ref. |44|). 
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The first moments of the Sivers and Boer-Mulders functions, i.e. the quantities 

<7 (1) (*) = jd 2 k T ^- 2 q(x,k T ) , (24) 

where q = hi Q , or f^p, are depicted in Figs. 2-5. The results for the calculation in the MIT 
bag model as well as in the CQM are represented for both functions on Fig. 2. As already 
mentioned in Refs. j^, EJ, the signs of the Sivers functions for the u and d flavors, nega- 
tive and positive respectively, are in agreement with both the theoretical and experimental 
knowledge. Concerning the Boer-Mulders function, are available but following the 

pattern predicted by quark helicity-flip IPD GPDs 30|, |31| , the same sign for the u and d 



flavors is expected. This is found in the present approach, confirming; the results of previous 
estimates, such as the first MIT bag model calculation 33J, the quark-diquark model 28j], a 
large N c analysis 34j and a phenomenological parameterization based on the quark-diquark 
picture 35|. This can be clearly seen in the lower panels of Fig. 2. The good performances 
of the approach, established in the calculation of the Sivers function, are therefore confirmed 
here for the Boer-Mulders function. 

For the sake of clarity, the same results are rearranged in Fig. 3, such that one can 
compare, in both models, the Sivers with the Boer-Mulders functions. It turns out that the 
shape of both functions is similar in the CQM and in the MIT bag model. In particular, 
the absolute value of /if turns out to be a little bigger than that of f^. Also the size 
is not too different, since it crucially depends on the scale of the model, which fixes the 
strong coupling constant g appearing in Eqs. ([TBI and (122]) . A small variation of g affects 
dramatically the results. We have decided to fix the hadronic scale as explained at the 
beginning of this section, i.e. in a model independent way, obtaining a s (/io)/(47r) ~ 0.13. 
One could have instead taken the initial scales corresponding to the ones fixed for the two 
different models by the Authors proposing them, i.e. a s {^) / '(47r) ~ 0.11 for the CQM 40] 
and a s (/ig)/(47r) ~ 0.18 for the MIT bag 43], respectively. The sizes of the results in the 
two different models would have come out more similar. Since there are no precise data 
available, a finer tuning of these parameters is not relevant at this moment. 

In Figs. 4 and 5, a microscopic analysis of the contributions of different helicity-flips of 
the quark interacting with the virtual photon and of the active recoiling one is presented, 
for fyp and respectively. 

In the case of the Sivers function, the helicity-flip can occur either for the quark interacting 
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with the virtual photon, or for the active recoiling one [4j, [5[ • The latter contribution has 



been disregarded in Ref. 



BQ. 



331 ]. while it is seen in Fig. 4 that it is very important. We recall 



that, once such a term is taken properly into account, the Burkardt Sum Rule is fulfilled to 
a large extent [5|. A similar problem occurs also for the calculation of h^. In the previous 
section, we have shown that, in the case of this function, contributions are found if there is 
no helicity-flip for the two active quarks, or if they both flip their helicities (cf. Eq. ( |T3|) 
and (120]) ). The latter contribution has been once more disregarded in Ref. 33|, while it is 
seen in Fig. 5 that it is sizeable. Moreover, if the two contributions are disregraded, f£p for 
a given flavor turns out to be proportional to hy for the same flavor. Besides, fyp (hi) for 
a flavor turns out to be proportional to (h^) for the other flavor. 

These proportionalities, are not found once the proper helicity-flip contributions are taken 
into account. Indeed, the reincorporated term is subtractive in the case of the Sivers function, 
while it is additive in the case of the Boer-Mulders function. This implies that the two T-odd 



functions cannot be proportional. In any case, the proportionality found in Ref. [33| has no 
physical motivation, as, in order to calculate the two T-odd functions, one is using a two- 
body operator associated with FSI and therefore one should not expect a proportionality 
between the u and d results. On the contrary, in the calculation of conventional PDs, in any 
SU(6) model calculation, the used operators are of one-body type and therefore the results 



turn out to be proportional 



39|]. 



Another interesting byproduct of the present results is that there is basically no quali- 
tative difference between the results in the two models, despite the fact that one of them 
is NR. As already said in the previous section, this has to do with the fact that the terms 
which encode the most relativistic contribution arising in the MIT bag model turn out to 
be a few orders of magnitude smaller than the dominant ones. This is seen in Eqs. (A8) 
and (A14) in the Appendix at the formal level. This gives us more confidence in the order 
0(k 2 /m 2 ) of the NR expansion used. 

Let us now compare our results with the other microscopic calculation of the Boer-Mulders 
function, the one in the quark-diquark model 281 ] . We get the same, negative, sign for the 
u and d flavors, which has been proven to be related to the model independent features of 
IPD GPDS 1 30]. On the other hand, the relative size of the u and d flavors turn out to be 

we have the result for u a bit larger than for d). Anyway, 



different (at variance with 



as far as we understand, a model parameter in Ref. 28J can change such a relative size, 
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making this comparison not conclusive. What makes our approach unique, transparent and 
instructive, is that a microscopic analysis of the helicity-flip of the different quarks involved 
is not possible in two-body models, such as the quark-diquark model. We have seen instead 
that one can learn several interesting aspects of the problem by performing this helicity-flip 
analysis at the single quark level. 

One relevant missing part in our analysis is the evaluation of the pQCD evolution of 
the results from the scale of the model to the experimental one. Unfortunately, to our 
knowledge, the anomalous dimensions of the relevant operators have not been calculated yet 



and no workable evolution equations have 



jeen provided. Some crucial steps towards the 



45 



46 



47 



48 



49 



50J, and they will be considered 



solution of this problem have been done 
in future work. 

Other natural developments of this calculation are the analysis in a relativistic three-body 
model, as well as the connection of T-odd TMDs with the corresponding IPD GPDs. 



IV. CONCLUSIONS 



We have applied a formalism, developed previously to calculate the Sivers function, to 
the calculation of the Boer-Mulders function with the aim of understanding the microscopic 
mechanisms associated with the transverse polarization of the nucleon. The interest of the 
calculation lies in the fact that, despite the difficult experimental extraction of this function, 
ongoing experiments will produce data. We have applied the formalism to two very different 
models of nucleon structure: a non relativistic constituent quark model and a relativistic 
MIT bag model, both of which are able to interpret well the static properties of the proton. 
Our calculation is complete in the sense that it takes into account all possible contributions, 
i.e. not only the helicity non-flip contribution, but also the helicity double flip one. 

The calculated moment of the Boer-Mulders function has the same negative sign for the u 
and d flavors confirming the results of previous predictions and estimates. This nice outcome 
also happened in our previous calculation of the Sivers function using the same formalism, 
where our results satisfied the Burkardt Sum Rule and followed the general trend of previous 
correct estimates. The shapes of both functions, Sivers and Boer-Mulders, are similar in the 
CQM and the MIT bag model. This is a consequence of the fact that the terms carrying the 
most relativistc contribution in the bag model are a few orders of magnitude smaller thean 
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the main contributions. It should be stressed that, the additional terms found in both of 
these functions in the MIT bag model calculations, are crucial to establish this behavior and 
also to satisfy fundamental relations like the Burkardt Sum Rule. 

To conclude, we have developed a formalism that allows a straightforward evaluation of 
all kinds of parton distributions, which can be easily applied to different models of hadron 
structure and, in particular, to two of the most succesful ones. Moreover the results ob- 
tained are quite model independent in structure, and are in agreement with fundamental 
principles and predictions. Our analysis has the advantage, with respect to other analyses 
in other formalisms, that it shows the microscopic separation of the different helicity-flip 
contributions for the different quarks involved in the process and, in so doing, is able to 
understand the microscopic structure of the transverse polarization up to evolution effects. 



APPENDIX A: SOME DETAILS OF THE RESULTS OF THE CALCULATION 

Some details are listed below. 

1. MIT bag model calculation 

The matrix elements Eq. (I12jl are listed below, for Q = u and S z = 1(— 1), 



C 



<++,++ 

liSz=l(-l) 




(Al) 




(A2) 




C 



(-+,-+ 

uS z =l(-l) 



0(0) 



(A3) 



and for Q = d and S z = 1(— 1), 



C 



<++,++ 

dS z =l(-l) 




(A4) 




a 



,+-,+- 

dS z = l(-l) 



0(0) 




(A5) 
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c ^s%U) = 1 (§) > c Z=iU) = (°) • ( A6 ) 

The functions entering Eq. ffT3l) are listed below: 
+ *f>*f>}\ 3 CAT) 



(A£ 

u f jsr /iw*i( fc 3). /,/x jxti(^), / Ml /- r V *i(^3) *i(fca) 1 
H 2 = - J d 6 k 3 Ucl -^—toih) ~ (*s - ? w ) -jr^to(h) + (^t x /c 3 J — ^ j . 



(A10) 



and 



7i = k J_M to{k/) _ (r _ ^ ^ to(fc ) _ C jfcV _ fc V) M^l M^l (All) 

Ar \ / k k' 

h = ky t -^t (k') + (ky- q y) t -^lt (k)+ ftV + JtV* 

k k' V > 



(A12) 



+ (tf* _ gxfcx _ fc «2 + qVtf + A;* 2 ) * l(A/) , (A13) 
v ' k k' 

J 2 = - (2k x k y - q x k y - q y k x ) tl( £ ' , (A14) 

with fc = and k! = k — (fa. 

The constant C entering Eq. (|T3j) is: 



9 (2tt) 

with 



C=\llSrz, (A15) 



^ = 4,^ = ^^ (A16) 



determined by the MIT bag model parameters 



1(3 



2. CQM calculation 



The potentials corresponding to gluon exchanges without quark helicity flips or with 
double quark helicity flips, Eq. f[2"Ul) . are listed below: 



yno-flip ft t ff\ . - _}_ 




q x q x k%\ 



2m ^ Am 2 J 



m 2m Am 2 J 3 



4m 2 J 



u 



3 I ^ | 1 \ ^zz 

1" o 1" i — 5" °3 a i 

m 2m Am 1 I 



(A17) 



\2m Am? J \ m 



i - <f fc 3 

4m 2 



+ H£ + SViH-s + ^y!S)^ n (A18) 



(fkf\ ( 
2m ' 4m 2 J \ m Am 2 

The set of conjugated intrinsic coordinates used in Eq. (l2"Tj) is 

R = —{fi + rl + r 3 ) ^ K = — (£ + fc 2 + fc 3 ) , 

A = -^(rl + rg - 2rl) <-> fc A = -^(fc* 2 + k 3 - 2k\) . (A19) 



The -B and A functions appearing in Eq. ( 12211 are 

B = -^L(l+ k P 3 ^ Q 2 x + <ly ~ QxK - q y k y 



2m \ 2 a/2 \/6 2m 4m 2 
A = - jj + + ^ (4g 3 - + 4^ - g^ 2 - q x q y k y - 3k x q 2 y ) (A20) 
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FIGURE CAPTIONS 



Fig. 1: The contributions to the Sivers and Boer-Mulders functions in the present approach 



The graph has been drawn using JaxoDraw 51]. 



Fig. 2: Upper (lower) panel, left: the first moment of the Sivers (Boer-Mulders) function, 
Eq. (12~4"1) . for the u flavor, in the NR (full) and MIT bag models (dashed); right: the same 
for the d flavor. 

Fig. 3: Left panel: the first moment for the Sivers and Boer-Mulders functions in the NR 
model. Dashed curve: the Sivers function for d; dot-dashed curve: the Sivers function for u; 
long dashed curve: the Boer-Mulders function for d; full curve: the Boer-Mulders function 
for u. Right panel: the same in the MIT Bag model. 

Fig. 4: Comparison between the results for the first moment of the Sivers function in the 
NR model and in the MIT bag model. Left panel: results in the NR model. Full (Long 
dashed) curve: full result for the flavor u (d); dot dashed (dashed) curve: result with the 



spin-flip occurring only for the interacting u (d) quark (as in Ref. 33[). Right panel: the 
same for the MIT Bag model 

Fig. 5: The same of Fig. 4, but for the Boer-Mulders function. 
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